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Abstract. Semiconductor superlattices are interesting for two distinct reasons: the 
possibility to design their structure (band-width(s), doping, etc.) gives access to a 
large parameter space where different physical phenomena can be explored. Secondly, 
many important device applications have been proposed, and then subsequently suc- 
cessfully fabricated. A number of theoretical approaches has been used to describe their 
current-voltage characteristics, such as miniband conduction, Wannier-Stark hopping, 
and sequential tunneling. The choice of a transport model has often been dictated 
by pragmatic considerations without paying much attention to the strict domains of 
validity of the chosen model. In the first part of this paper we review recent efforts 
to map out these boundaries, using a first-principles quantum transport theory, which 
encompasses the standard models as special cases. In the second part, focusing in the 
mini-band regime, we analyze a superlattice device as an element in an electric circuit, 
and show that its performance as a TIfz-photon detector allows significant optimiza- 
tion, with respect to geometric and parasitic effects, and detection frequency. The key 
physical mechanism enhancing the responsivity is the excitation of hybrid Bloch-plasma 
oscillations. 



1 Introduction 

Ever since the pioneering work of Esaki and Tsu [1], which drew attention to 
the rich physics and potential device applications of semiconductor superlattices, 
these man-made structures have remained a topic of intense research. Semicon- 
ductor superlattices have proven to be a fruitful platform for studying a wide 
range of transport phenomena, such as their intrinsic negative differential con- 
ductivity [2] , the formation of electric field domains [3] , Bloch oscillations [4] , as 
well as dynamical localization [5] and absolute negative conductance [6] under 
external irradiation, just to mention a few. 

These phenomena depend crucially on the relations of the energy scales in- 
volved, namely the zero-field miniband width (which is four times the interwell 
coupling Ti), the scattering rate F/h^ and the potential drop per period (= eFd^ 
where F is the applied static field and d is the superlattice period). Three distinct 
approaches have been used to describe transport in the parameter space spanned 
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by {Ti,eFd,r): miniband conduction (MBC)[1], [7], Wannier-Stark hopping 
(WSH)[8], and sequential tunneling (ST) [9], [10]. Until recently, however, the 
precise range of validity of these various approaches had not been addressed 
quantitatively. To achieve this, one much use a first-principles approach, which 
reduces to the standard theories in the appropriate limits. After a brief review of 
the standard approaches in Sect. II, we introduce in Sect. Ill a nonequilibrium 
Green function formalism, which we have used to delineate the boundaries of 
the different domains of validity [12], [13]. We also present comparisons with 
Monte Carlo simulations. As we shall see, under favorable conditions it is quite 
possible to obtain a very accurate description with the standard methods, which 
is a great advantage because the first-principle Green function calculations are 
quite complicated and have so far successfully implemented only for rather sim- 
ple scattciring nicichanisms (the scattering matrix elements arc assumed to be 
independent of momentum transfer). Thus, in Sect. IV we adopt one of the 
standard approaches, i.e., the miniband approach, to model a superlattice THz- 
photon detector, taking into account the effects due to the external circuitry. We 
conclude that by detailed modeling substantial device performance optimization 
can be achieved, e.g. the detector sensivity may be improved by almost 50 % by 
a judicious choice of its parameters. 

2 Standard transport models 

Here we review the standard models used to describe transport in semiconductor 
superlattices. For simplicity, the results quoted in the next three subsections arc 
written in a relaxation time approximation, but a generalization to more realistic 
scattering processes is possible. As an underlying Hamiltonian we use a tight- 
binding model: 

Hn^m = iKm-1 + ^n,m+l) + ^„,„(^fc - ncFd) . (1) 

Here Ti is the overlap matrix element, E/. = h^k'^/{2m) is the kinetic energy 
perpendicular to the growth direction, and the electric field is taken into account 
by a shift in the site energies. 

2.1 Miniband conduction (MBC) 

For zero electric field Eq. (1) is diagonalized by a set of Bloch fimctions ipq{z) = 
gin9d^^^_-,-j (^iigpg ^iig wave function if'nl-z) is localized in quantum well n, 
e.g., a Wannier function) and the dispersion relation is given by the miniband 
E(q) = 2TiCos{qd). The stationary Boltzmann equation for the distribution 
function f{q, k) is then 



eF dfjq, k) _ nFjEjq) + E^) ~ f{q, k) 
h dq T{E{q) + Ek) 



(2) 
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where the relaxation-time is allowed to depend on energy, see, e.g., Ref.[12] for 
a suitable model. Once the solution to Eq. (2) is found, the current is calculated 
from 

The electron density per period is given by 



N2D = ^ f d^k £^\f{q,)<.) (4) 



and is used to determine the chemical potential for a given electron density. This 

approach can be extended beyond the relaxation time approximation[14], [15], 
but the generic features remain unchanged. 



2.2 Wannier-Stark hopping (WSH) 



In the presence of an electric field, the eigenstates of the Hamiltonian become 
the localized Wannier-Stark states. 



rt ^ 



2Ti 



(5) 



with energy E^, = —veFd, and is the Bessel function of the first kind. 

Scattering causes hopping between the different states. Within Fermi's golden 
rule, the current is given by 



(6) 



Here the term YJ\.'^'nJn-i^ arises due to the spatial overlap of the Wannier-Stark 
functions. The electron density per period is given by: 



= PoksTe log 



(7) 



which relates /v, to N2d- Again, it is possible to generalize this approach to more 
realistic scattering mechanisms [16], [17]. 
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Fig. 1. Overview of standard theoretical models for superlattice transport 



2.3 Sequential tunneling (ST) 

In this approximation the phase information is lost after each tunneling event 
between adjacent wells. The scattering within a well is treated self-consistently by 
solving for the spectral functions A{£, k); in this work we use the self-consistent 
Born-approximation[18] for the self-energy. The transitions to neighboring wells 
are calculated in lowest order of the coupling yielding [20], [18], [19]: 

J{F) = ^ [ d?k f ^T^Ai£,k)A{£ + eFd,k)[nF{£)-nFi£ + eFd)] . 

ZTT J J ZTTft 

(8) 

The carrier density is given by: 

^ i/ "^'^ / ^"^'(^)^(^'k) . (9) 

This approach gives quantitative agreement with experiments in weakly cou- 
pled structures when realistic models for impurity and interface scattering are 
employed [10], [20]. 

2.4 Summary 

The important issue to recognize is that these three approaches treat scatter- 
ing, external field, and coupling within different approximations. MBC does not 
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properly include field-induced localization because of its inherent assumption of 
extended states, WSH treats scattering in lowest order perturbation theory (in 
particular, there is no broadening of the states), and ST is explicitly lowest order 
in the interwell coupling. The basic features of these models are summarized in 
Figure 1. 

2.5 Nonequilibrium Green functions (NGF) 

The basic building blocks arc the correlation and retarded Green functions: 

G<Jt,f',k) = «(at(t',k)a™(t,k)) (10) 
Glfjt,t',k) = -te{t-t'){{a^{t,-k),ai{t',k)}), (11) 

where and al^{t,k) and a„(t,k) are the creation and annihilation operators for 
the state 'Fn{z)e^^^''^y A in well n. These functions obey the Dyson and Keldysh 
equations, respectively, given below for the superlattice case. Ref.[21] may be 
consulted for a text-book discussion. The observables, such as the momentum 
distribution, current, and electron density are computed from ^{t,t',k): 

/QM(g, k) = / ^ E ''~"""'G<,{£, k) (12) 

n 

^(^^ = d'klm.{TiG<+i,n{^,k)} (13) 

These expressions exploit the fact that in the stationary state the Green functions 
only depend on the time difference r = t — t', and one can define a Fourier 
transformation via [21]: 

G„,„(f,k) = JdTe<^-'^'''^)^GmAt,t-T,k). (15) 

Without scattering between the k-states and at Ti = the Green-functions are 
diagonal in the well index: G^*„(£,k) = 5m,n9^^{^,^) with the free particle 
Grccn-function ,gjf*(f ,k) = l/{£ — Ek + iO+). The full Green function is then 
determined by the Dyson equation: 

{£,k)= gZ' {£, k) + gZ' + eFd^^^k^ 
xY^^Z + efd^,k) Gl% ^£ + eFd^-^,k^ . (16) 
The self-energy will be written as 

^m%i^^ ^) = ^m,n^™'(f, k) -|- Ti6m+l,n + Ti6m-l,n (17) 
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where S contains contributions both from impurity and phonon scattering. The 
relevant expressions are 

^n^^i^) = f E^-pG",n^<(^'k') (18) 
k' 

for impurities, while for optical phonon scattering we take (see, e.g., Ch. 4.3 of 
Ref. [21] for the derivation): 

4<o(f) = ^ E {^oG< „(£ - Tico, k') 

k' 

+ (7Vo + l)G<„(f + ;iu;o,k')} (19) 



SZ^£) = ^E{(^o + l)^tn(^-?^'^o,k')+iVoG-;„(£ + ?ic^o,k') 

k' 

It is possible to simulate acoustic phonons with a similar expression, using a 
small fictitious discrete energy TiuJac [13]. In the above, we also gave the self- 
energy expressions needed in the Keldysh equation: 

G< ,„(f , k) = ^ GZ]m. U + eFd"^^, k) G^t„ (s + eFd^^^, k) 

xE<f£ + eFd(m,-'^),k) . (21) 



The numerical evaluation of these equations has been discussed in two recent 
publications [12], [13], and here we just summarize our basic strategy, and give 
a few representative results. The first step in the analysis consists of evaluating 
the current-voltage characteristics for the different models, and an example of 
such a calculation is given in Fig. 2. We note that the Boltzmann equation gives 
an excellent agreement with the full quantum result, in particular for low electric 
fields. The additional structure in the quantum mechanical curve is a phonon- 
replica: the Boltzmann equation cannot capture features like this because the 
electric field docs not enter as an energy-scale in the collision integral. It is also of 
interest to note that the use of a realistic electron-phonon scattering model, in- 
stead of a simple relaxation time, leads to a cur rent- volt ago characteristic which 
deviates significantly from the simple Esaki-Tsu shape, Wdrift oc F/(F^ + F^^i-jJ. 
By performing a large number of calculations like the one depicted in Figure 2 
it is possible to develop visual criteria as to the validity of the various simpli- 
fied approaches discussed in Section 2. These phenomenological considerations 
can be made more rigorous by examining in detail the retarded Green function, 
which determines the quantum mechanical correlation between quantum wells 
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Fig. 2. Drift velocity for a wide-band superlattice. Pull line: Calculation by nonequilib- 
rium Green functions. Dashed line: Monte Carlo simulation of Boltzmann's transport 
equation (from Ref.[12]). 



n and m [12]. For example, if the terms G53:i n the order of G^n\n the 

wave- function is dclocalizcd. By using the expressions given in Ref. [12], one finds 
that the states are delocahzed if 2|Ti| ^ F and 2|Ti| ^ eFd, i.e., the Boltzmann 
miniband picture is a useful starting points. Similarly, by demanding that the 
^Inn vanish if m 7^ n, one can find the regime where sequential tunneUng dom- 
inates. The respective ranges of validity can be presented as a "phase digram" , 



Additional insight to the differences between quantum kinetics and Boltz- 
mann kinetics can be obtained by examining the momentum distribution func- 
tions. Figure 4 presents such a comparison. We direct attention to the following 
features. At low fields. Fig. 4(a), the distribution function can be viewed as a 
distorted thermal equilibrium function, and linear response theory holds, cor- 
responding to the linear part of the current-voltage characteristics of Fig. 2 
(however, at very low temperatures one should consider weak localization ef- 
fects, not included in the present choice of the impurity self-energy, and the 
agreement between quantum and Boltzmann calculations may be weakened). If 
the electric field increases, electron heating becomes important. For moderate 
fields the distribution function resembles a distorted equilibrium function, but 
with an elevated electron temperature. Fig. 4(b). At even higher fields, close 
to the maximum in the IV-curve, the distribution function strongly deviates 
from any kind of equilibrium in q space. Fig. 2(c). The results for quantum and 
semiclassical calculations look similar in the negative differential conductance 
(NDC) regime (not shown in Fig. 4.). The situation changes dramatically when 
the field energy is larger than the mini-band width, edF > A, Fig. 4(d). As the 
electrons can perform several Bloch oscillations in the semiclassical picture, the 
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Fig. 3. Ranges of validity for different transport models 



distribution function is almost flat withing the Brillouin zone of the miniband. 
The latter holds for the NGF result as well. However, the absolute values of the 
distribution functions differ significantly. The reason is the modification in scat- 
tering processes due to the presence of the electric field, leading to significant 
deviations in the distribution function. We noticed this phenomenon already in 
our discussion of the current-voltage characteristics of Fig. 2. The strong changes 
in the quantum momentum distribution are further illustrated in Fig. 5. 

We can summarize the results given in the first part of this paper as follows. In 
wide-band superlattices the Boltzmann equation gives reliable results concerning 
linear response at low fields, electron heating at moderate fields, and the onset 
of negative differential conductivity. In contrast, for high electric fields or weakly 
coupled SLs significant differences appear. In this case the quantum nature of 
transport is important and a semiclassical calculation may be seriously in error. 

3 Superlattice as a THz-photon detector 
3.1 Introduction 

It has recently been estimated [22] that the room temperature current responsiv- 
ity of a superlattice detector ideally coupled to the THz-photons can nearly reach 
the quantum efficiency e/fnj in the limit of high frequencies uj ^ v (here m is the 
incident radiation frequency and v is a, characteristic scattering frequency) . This 
value of the responsivity is being normally considered as a quantum limit for 



Transport in semiconductor superlattices 9 




Fig. 4. Electron distribution versus quasirnorncnturn in the Brillouin zone of the mini- 
band for different values of k (from Ref.[12]). The parameters as Fig. 2. Full line: NGF 
calculation. Dashed line: MC-simulation of BTE. The dotted line shows the thermal 
distribution oc exp[-S(g, k)/jfcsTe] with Te = T in (a) and Te = 140 K in (b) for 
comparison. 




Fig. 5. Electron distribution versus quasimomentum k (from Ref.[12]). (a,b,c): Com- 
parison between NGF calculation (full line) and BTE (dashed line); (d): Results from 
NGF for different fields. The BTE result (not shown) resembles the result from (c) for 
all three fields. 
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detectors based on superconducting tunnel junctions operating at low tempera- 
tures [23] . For high frequencies the mechanism of the THz-photons detection in 
superlattices was described [22] as a bulk superlattice effect caused by dynamical 
localisation of electrons. 

Here we describe a recent theory of the superlattice current responsivity [24] . 
We focus on relative broad-band superlattices {A ~ 20 meV), and field strengths 
smaller than the onset of negative differential conductivity. Thus, according to 
the analysis given above, a Boltzmann equation based description should suffice. 
One should note, however, that here we investigate a time-dependent situation, 
and that a microscopic analysis in the spirit of the first part of this paper has 
not, to our knowledge, been carried out. We use an equivalent circuit for the 
superlattice coupled to a broadband antenna (see Fig. 6), which is similar to 
the equivalent circuit used in resonant tunneling [25] and Schottky diode [26] 
simulations. The suggested equivalent circuit of the device allows one to treat 
microscopically the high-frequency response of the miniband electrons and, si- 
multaneously, take into account a finite matching efficiency between the detector 
antenna and the superlattice in the presence of parasitic losses. Our analytic re- 
sults lead to the identification of an important physical concept: the excitation 
of hybrid plasma-Bloch [27] oscillations in the region of positive differential con- 
ductance of the superlattice. Numerical computations [24] , performed for room 
temperature behavior of currently available superlattice diodes, show that both 
the magnitudes and the roll-off frequencies of the responsivity are strongly influ- 
enced by this effect. The excitation of the plasma-Bloch oscillations gives rise to a 
resonant-like dependence of the responsivity on the incident radiation frequency, 
improving essentially the coupling of the superlattice to the detector antenna. 
We will also show that peak current densities in the device and its geometrical 
dimensions should be properly optimized in order to get maximum responsivity 
for each frequency of the incident photons. Finally, we will present numerical 
estimates of the responsivity for the 1-4 THz frequency band and compare its 
value with the quantum eflaciency e/hw of an ideal detector. 

3.2 Theoretical formalism 

The exact solution of the time-dependent Boltzmann equation in the relaxation 
time approximation for an arbitrary time-dependent electric field can be pre- 
sented in the form of a path integral [28] : 



Using Eqs.(3) and (22) we find the time-dependent current I{t) describing ac 

transport in a superlattice. with electron performing ballistic motion in a mini- 
band according to the acceleration theorem and suffering scattering [29] , [30] : 




(22) 
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Fig. 6. (a) THz-radiation coupled to a A'^-pcriod semiconductor superlattice by a 
co-planar broad band bow- tie antenna, Pi and Pr arc the incident and reflected powers 
respectively, (b) Equivalent circuit for a THz-photon detector with a dc voltage bias 
source: S-miniband electrons capable to perform Bloch oscillations, C-superlattice car 
pacitance, iZs-parasitic series resistance, ^A-bow-tie antenna impedance, Vdc-dc bias 
voltage. 



where y(t) = LF{t) is the voltage across the superlattice perpendicular to the 
layers, L = Nd is the superlattice length, N is the number of periods in the 
superlattice sample, Ip = Sjp, S = wa^ is the area of the superlattice, a is the 
superlattice mesa radius, and 



vq f 2dkdq , ,s . , , x 



(24) 
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is the characteristic current density. The integration over q in Eq.(24) must be 
carried out over the Brillouin zone —-K/d < q < n/d. 

The peak current density jp and the scattering frequency v can be considered 
as the main parameters of the employed model. They can readily be estimated 
from experimentally measured or numerically simulated values of Ip and Vp. 
For both degenerate and non-degenerate electron gas one gets [29] , [30] 

Jp = en— (25) 

if A ^ kT,€F, where kT is the equilibrium thermal excitation energy, ep = 

?i^(37r^iV2£))^/''/(2mofr) is the Fermi energy of degenerate electrons, rricn = 
rnjzvn^l'^ is the density of states effective mass near the miniband bottom, and 
rrizz = 2h'^/Ad'^ is the effective mass of electrons along the superlattice axis. In 
the particular case of the Boltzmann equilibrium distribution function Eq.(24) 
yields [30] jp = {envo/2)[Ii{A/2kT)/Io{A/2kT)], where Io,i are the modified 
Bessel functions. 

Wc now suppose that in addition to the dc voltage Vsl, an alternating sinu- 
soidal voltage with a complex amplitude is applied to the superlattice: 

V{t) =VsL + l K exp(ia;i) + exp(-ia;i)] . (26) 

Generally, Vsi, K; can be found from an analysis of the equivalent circuit given 
in Fig. 6. We write the ac voltage amplitude as 14, = |K,|e''''; both \Vu,\ and ip 
can be obtained self-consistcntly taking account of reflection of the THz photons 
from the superlattice and their absorption in the series resistor Rs- 
Making use Eq.(23) we obtain [30]: 

/■oo 

I{t) = 2Ip I i/rffi exp(— i/ti) sin 

where 

^{t,h) = X ^ {iK,exp(ia;i)[exp(-za;ti) - 1] + c.c.} . (28) 

According to Eq.(27), electrons in a superlattice miniband perform damped 
Bloch oscillations with the frequency = eVsh/Nh = eEshd/h, and the 
phase ${t,t\) modulated by the external ac voltage. 

Equation (27) contains, as special cases the following results: (i) a harmonic 
voltage V{t) {Vsl = 0) leads to dynamical localisation, and current harmonics 
generation with oscillating power dependence [30] ; (ii) a dc current- voltage char- 
acteristics of the irradiated superlattice lociYsLjVu)) = (w/27r) / I{t)dt shows 
resonance features ('Shapiro steps') leading to absolute negative conductance 
[30], [31], [32], [33]; (iii) and to generation of dc voltages (per one superlattice 
period) that are multiples of /ia)/e[34]. 



eVsL 



(27) 
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3.3 Current responsivity 

We define the current responsivity[23] of the superlattice detector as the current 
change AI induced in the external dc circuit per incoming ac signal power Pj: 

Ri{co,VsL) = ^ . (29) 

This definition takes into account both the parasitic losses in the detector and 
the finite efficiency for impedance-matching of the incoming signal into the su- 
perlattice diode. 

It can be shown [24], that in the small-signal approximation both the dc 
current change AI^^ and the power Pj^^ absorbed in the superlattice are pro- 
portional to the square modulus of the complex voltage |VL p. This circumstance 
permits us to calculate |K;p self-consistently for given values of the incoming 
power, making use a linear ac equivalent circuit analysis and, then, find the 
current responsivity Ri{Lo, Vsl)- 

The results of the calculation of the superlattice current responsivity Vsl) 
are presented in the following form: 

Ti(nV ^_ Rf\u^,VsL)A{u:,VsL) 

"^^^ - l + Rs{dll,UVsL)/dVsL) ' ^'"^ 

where 

Ti^^)(nV \ - ^ {ysL/Vp)[^ + [u^rf - [Vsl/Vp?] 

K^,vsL) Nnv[l + {VsL/VpY][l + {ojTY-{VsL/VpY] ^ ' 

is the superlattice current responsivity under conditions of a perfect matching 
and neglecting parasitic losses, Rs 0. [22] 

The factor A{lu, Vsl) in Eq. (30) describes the effect of the electrodynamical 
mismatch between the antenna and the superlattice and the signal absorption 
in the series resistance 



A{uj,Vsl) = 



ZA-{Zi^{co,VsL) + Rs) 



ZA + {Zi^{u,VsL) + Rs) 



RcZi^,{u;,VsL) 



ReZ^L^iv,VsL) + Rs ' 

(32) 

The first factor in Eq. (32) describes the reflection of the THz-photons due to 
mismatch of the antenna impedance Za and the total impedance of the device 
Z^q{lj,Vsl) + Rs, with the second one being responsible for sharing of the 
absorbed power between the active part of the device described by the impedance 
^Aci^^ ^sl) and the series resistance Rs- 
The superlattice impedance is defined as 

Zi^ito, Vsl) - 1/ [G^Aci^. Vsl) + tcoC] , (33) 

where G^q{u), Vsl) is the superlattice conductance, C = cqS/A-kL is the capac- 
itance of the superlattice, and eo is the average dielectric lattice constant. 
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Finally, the last factor in the denominator of Eq. (30) describes the redistri- 
bution of the external bias voltage Vdc between the dc differential resistance 
of the superlattice {dI^^{VsL) / dVsL)~^ and the series resistance Rs, with the 
dc voltage drop on the superlattice Vsl being determined by the solution of the 
well-known load equation [23] 

Vdc = Vsl + /^c {Vsl)Rs ■ (34) 



3.4 Superlattice dielectric function. Hybridisation of Bloch and 
plasma oscillations 

We next analyze the condition of optimized matching of the superlattice to the 
incident radiation. Assuming the limit of negligible series resistance Rs this 
condition can be obtained from the solution of the equation 

Z^^iiJ,VsL) = Za (35) 

for the complex frequency u){Vsl) ■ This solution determines the resonant line 
position and the line width at which the absorption in the superlattice tends to 
its maximum value. 

One can transform Eq. (35) to the following form: 

<'^'^-) = ^' ^''^ 
where ^ 

e{Lo,EsL) = £0 + -^FiiLO,EsL) (37) 

is the dielectric function of the superlattice, with the dc field Esl being applied 
to the device [26], and Fi{u),Esl) is defined by 

F(u.V ) l + icOT-jVsL/Vp)^ 

' [1 + {Vsl/Vp)^] [(1 + icoT^ + {VsL/Vpn ■ ^ ' 

In the high-frequency limit €q/CZaco — *■ the solution of Eq. (35) coincides 
with the solution of the equation 

eiLu,EsL) = (39) 

describing the eigenfrequencies uj^ of the hybrid plasma-Bloch oscillations in a 
superlattice, [26] 

-,1/2 



w±(£'sl) = ±wp 



1 / 



1 + {Esl/Ep)^ \wpj 



+ — [Esl/EpY 



+ iu (40) 



where wp is the plasma frequency of electrons in a superlattice. The plasma 

frequency ujp can be given in terms of the small-field dc conductivity ctq or, 
equivalently, in terms of the peak current density jp 
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Fig. 7. The calculated hybrid plasma-Blocli oscillation frequency fn as a function of 
the normalized superlattice voltage drop Vsl/Vp for different values of the peak current 
densities jp = 10, 50, 100, 300, 500, and 1000 kA/cm^ (from Ref.[24]). Typical values 
of the superlattice parameters { d = 50A, Ep=W kV/cm, eo =13) were used for the 
calculations. 



Equation (41) reduces in the particular case of wide-miniband superlattices 
{A ^ kT^ep) to the standard formula ojp = {^-Ke^n/ eomzzY^'^ ■ 

In the limiting case of small applied dc electric fields Esl I Ep one finds 
from Eq. (40) the plasma frequency — > ±u}p, while in the opposite case 
Egh/Ep — > 00, the Bloch frequency — > ±i7s = ±eEsLd/fi is recovered. 
The scattering frequency v in Eq. (40) is responsible for the line width of the 
plasma-Bloch resonance. 

We have calculated the hybrid plasma-Bloch oscillation frequency fn = 
0)^/2^, using Eqs. (40) and (41), for the typical values of the superlattice pa- 
rameters eo ~ 13, d ~ 50 A, Ep 10 kV/cm, = u/2n = 1.2 THz for different 
values of the current densities jp (see Fig. 7). For small values of the current 
densities jp ~ 10 kA/cm^ the frequency of the hybrid oscillation inc;reases with 
applied voltage in all range of the parameter Vsl/Vp . On the other hand, 
for higher values of the current densities jp ~ (50 — 1000) kA/cm^ the hybrid 
oscillation's frequency starts to decrease with increasing bias voltage in the sub- 
threshold voltage range Vsl < Vp. Then, at super-threshold voltages Vsl > Vp, 
u)H starts to increase again tending to the Bloch frequency. It is important to 
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note that at high values of the dc current densities jp the hybrid plasma-Bloch 
oscillations become well defined eigenmodes of the system {fn > fv)- There- 
fore, an essential improvement of the matching efficiency between antenna and 
the superlattice can be expected in the high- frequency range due to a resonant 
excitation of this eigenmode in the device. 
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Fig. 8. Real space energy diagram illustrating THz- photon (/ ^ v/2n) detection 
in the superlattice: DC electric field Esl is applied to the iV-period semiconductor 
superlattice with the miniband width A . Under the action of the dc field electrons 
perform Bloch oscillations with the spatial amplitude A/eEsL ■ At critical dc electric 
voltage (field) Vsl ~ Vp = NTiv/e {eEshd — hv) electrons move against the dc electric 
force due to absorption of photons climbing up the Wannier- Stark ladder. The energy 
2eVp should be absorbed from external ac field in order to subtract one electron from 
the external circuit. One half of this energy is needed for the electron to overcome the 
potential barrier which is formed by the dc force, with the other half being delivered 
to the lattice due to energy dissipation. A quasi-classical description of the process is 
valid if / <C A/U when allowed transitions between different Wannier-Stark state exist. 



3.5 High-frequency limit 

Let us compare the high-frequency limit of the responsivity of the superlattice 
with the quantum efficiency i?max = e/Tujj which is believed to be a fundamen- 
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tal restriction for the responsivity of superconductor tunnel junctions [23]. This 
quantum efficiency (or quantum hmit) corresponds to the tunnelhng of one elec- 
tron across the junction for each signal photon absorbed [23], with a positive sign 
of the responsivity. 

In our case tlic mechanism of the photon detection is different (see Fig. 
8). Electrons move against the applied dc electric force due to absorption of 
photons. At VsL = Vp the responsivity is negative, indicating that one electron 
is subtracted from the dc current flowing through the superlattice when the 
energy 2eVp is absorbed from the external ac field. One half of this energy 
is needed for the electron to overcome the potential barrier which is formed 
by the dc force, with another half being delivered to the lattice due to energy 
dissipation. If the applied dc voltage is strong enough, i.e. Vsl ^ Vp, dissipation 
plays no essential role in the superlattice responsivity. In this case the energy 
eVsL should be absorbed from the ac field in order to subtract one electron from 
the dc current simply due to the energy conservation law. 

3.6 Excitation of the plasma-Bloch oscillations 

For demonstration of the frequency dependence of the superlattice current re- 
sponsivity in the THz-frequency band we will focus on the GaAs/Gao.sAlo.sAs 
superlattices specially designed to operate as millimeter wave oscillators at room 
temperature. In Ref. [35] wide-miniband superlattice samples with d = 50 A, 
A 2± 113 meV, n ~ 10^^ cm~'^, were investigated experimentally. They demon- 
strated a well-pronounced Esaki- Tsu negative differential conductance for Esl > 
Ep ~ 4 kV/cm with the high peak current of the order of jp ~ 130 kA/cm^. 
The measured value of the peak current is in a good agreement with the esti- 
mate jp ~ (80 - 160) kA/cm2 for n ~ (1 - 2) x 10", T = 300 K based on Eq. 
(24), if one assumes an equilibrium Boltzmann distribution for the charge carri- 
ers. From the peak electric field and current we find the scattering and plasma 
frequencies fi, ~ 0.5 THz, fp = 2 THz, respectively, assuming eg = 13 for the 
average dielectric lattice constant. The maximum frequency for the semiclassical 
approach to be valid for these samples is Ja — 27 THz. 

Figure 9 shows the frequency dependence of the normalised current respon- 
sivity calculated for three values of the peak current density in the superlattice, 
i.e. jp = 13, 130, and 300 kA/cm^ and for three values of the peak electric field, 
Ep = 4, 9, and 13 kV/cm. We also use the typical values for the superlattice 
length L — 0.5/im (superlattice consists of 100 periods), and assume a ~ 2/im 
for the superlattice mesa radius [35], [36], [37]. We choose Rs = 10i7 for the 
series resistance of the device in the THz- frequency band, i.e. the same value 
as for resonant tunnelling diodes having the same radius of mesas [25]. The cal- 
culations are performed in the region of the positive differential conductance for 
dc bias voltage close to the peak voltage {Vsl = 0.95 Vp) 

For Ep = 4 kV/cm (/^ ~ 0.5 THz) Fig. 9 demonstrates well-pronounced 
resonant behavior of the normalised responsivity as a function of frequency. The 
resonance frequency and the maximum value of the responsivity rise if the peak 
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f(THz) 

Fig. 9. The frequency dependence of the normalised current responsivity 

|-Rijv| = \Ri/{c/T^^)\ of the supcrlatticc THz-photon detector (a = 2/im, L — 0.5/im, 
Rs = 10i7, VsL ~ 0.95Vp) for three values of the peak current density {jp = 13, 30, 
and 300 kA / cm'^ ) and for three values of the peak electric field {Ep — 4, 9, and 13 
kV/cm ) (from Ref.[24]). The relevant positions of the hybrid plasma-Bloch frequencies 
fn are indicated for each curve by arrows showing characteristic resonance (high peak 
current densities) and roll-off (low peak current densities) behavior. 



current density increases. For jp = 300 kA/cm^ the normalised responsivity 
reaches its maximum value —RiN — 0.02 (— i?i ~ 2 A/W) at frequency / ~ 2.5 
THz. For higher values of the peak electric fields Ep = 9 kV/cm {f^, ~ 1.08 
THz) and Ep = 13 kV/cm (/^ ~ 1.57 THz) the resonance line-widths are 
broadened due to implicit increase of the scattering frequencies. In particular, 
for Ep = 13 kV/cm, jp ~ 300 kA/cm^ the normalised responsivity has an 
almost constant value —RiN — 0.006 {-Ri ~ 0.6 A/W) up to / ~ 2.5 THz and, 
then, rapidly decreases. The frequency behavior of the normalised responsivity 
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originates from excitation of the plasma- Bloch oscillations in the superlattice. 
We indicate in Fig. 9 the positions of the hybrid frequencies fu = |/27r with 
arrows. For small peak electric fields (low values of the scattering frequencies) the 
hybrid frequency corresponds to the maximum of the normalised rcsponsivity. 
For higher values of the peak field (higher values of the scattering frequencies) it 
corresponds to the roll-off frequency at which the responsivity starts to decline. 



3.7 Optimized superlattice length 

The cinhanccmcnt of the normalised responsivity requires an optimum matching 
efficiency of the superlattice to the broad-band antenna and minimization of the 
parasitic losses in the series resistor. These requirements impose an optimum 
length of the superlattice for each chosen frequency of the incoming THz-photons 
and series resistance. 




Fig. 10. The dependence of the normalised current responsivity |-RiAr| = \Ri/{e/huj)\ 
of the superlattice THz-photon detector (a = 2/im, Ep = 4kV/cm, jp = 130 kA/cm^) 
at / = 2.5 THz for three values of the series resistance (-Rs=10, 30, and 50 Q) as a 
function of number of the superlattice periods N (from Ref.[24]). 



We show in Fig. 10 the dependence of the normalised responsivity on the 
number of the superlattice periods for / = 2.5 THz. We used for calculation 
a = 2/zm, jp = 130 kA/cm^, Vsl = 0.95 Vp, and three values of the series 
resistance Rs — f 0, 30 and 50 Q. For all three values of the series resistance the 
responsivity displays a well pronounced maximum for the optimum number of 
the superlattice periods A'' = iVmax- The value of Amax increases with increasing 
of the series resistance (ATmax — 40 for Rs = 10 i7, A^max — 60 for Rs = 30 
Q and for Amax — 90 for Rs = 50 Q). This result can be readily understood 
by recalling that a larger volume of the superlattice minimizes parasitic losses 
for higher values of the series resistance because of reduction of the sample's 
capacitance. 
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In addition to the examples discussed here, Ref. [24] reports several other as- 
pects of the superlattice responsivity, such as bias voltage dependence, optimized 
peak current density etc. 

3.8 Conclusions 

We have illustrated here the steps required to perform a superlattice device op- 
timization. These include: (i) Consideration of the frequency dependence of the 
superlattice. Here a truly microscopic analysis is not yet available, but a Boltz- 
mann equation based theory should be applicable, with due caution, (ii) The 
impedance matching between the antenna and the superlattice must be opti- 
mized, (iii) The effect of parasitic losses must be included in the analysis. One of 
the main results emerging from our analysis is the importance of collective exci- 
tations: the hybridized plasma-Bloch oscillations both enhance the responsivity 
and increase the roll-off frequency. It is found that an optimized superlattice can 
have a responsivity which approaches 10 % of an ideal superconducting tunnel 
junction. The great advantage, of course, is that superlattice based detectors 
work even at room temperature, and this property should offer a wide range of 
applications. 
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